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ABSTRACT 

To model the dynamics of the chiral order parameter in a far from equihbrium 
QCD phase transition, we consider quenching in the 0(4) hnear sigma modeL 
We summarize arguments and numerical evidence which show that in the period 
immediately following the quench arbitrarily long wavelength modes of the pion 
field are amplified. This results in large regions of coherent pion oscillations, and 
could lead to dramatic phenomenological consequences in ultra-relativistic heavy 
ion collisions. 

My talk was a description of work done with Frank Wilczek [1]. In these 
proceedings, I sketch our central results, emphasizing several points that were 
raised in discussions at the conference. The interested reader should, however, 
consult Ref. [1] and our earlier work [2] for a more detailed exposition. 
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1. Misalignment of the chiral condensate 



Among the most interesting speculations regarding ultra-high energy hadronic 
or heavy nucleus coUisions is the idea that regions of misahgncd vacuum might 
occur [3]. Misaligned regions are places where the four-component field 0" = (a, vf), 
that in the ground state takes the value {v, 0) is instead partially aligned in the 
7? directions. Because of the explicit chiral symmetry breaking {i.e because the 
pion is not massless), in such a region would oscillate about the a direction. 
If they were produced, misaligned vacuum regions would relax by coherent pion 
emission — they would produce clusters of pions bunched in rapidity with highly 
non-Gaussian charge distributions. In each such cluster, the ratio 



is fixed. Among different clusters, R varies and is distributed according to 



As an example of (1.2), we note that the probability that the neutral pion fraction it! 
is less than .01 is .1! This is a graphic illustration of how different (1.2) is from what 
one would expect if individual pions were emitted with no isospin correlations many 
pions. We have proposed [1] a concrete mechanism by which such phenomena may 
arise in heavy ion collisions for which the plasma is far from thermal equilibrium. 



R = 



n^o + + n. 




V{R) = \r~^''^ . 
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2. Emergence of long wavelength pion 
oscillations following a quench 



In studying the behaviour of the plasma in the central rapidity region of a heavy 
ion collision at RHIC energies or higher, it seems reasonable to assume that after 
a time of order 1 fm a hot plasma is formed in which the chiral order parameter is 
disordered and in which the baryon number density is low enough that it can be 
neglected. Our goal is to study the behaviour of the long wavelength modes of the 
chiral order parameter as this plasma loses energy and a develops an expectation 
value. 

In previous work (Rcf. [2], references therein, and Wilczck's talk at this con- 
ference) we considered the equilibrium phase structure of QCD. We argued that 
QCD with two massless quark flavours probably undergoes a second order tran- 
sition. At first sight, this might seem ideal for the development of large regions 
of misaligned vacuum, since the long wavelength critical fluctuations characteris- 
tic of a second order transition are such regions. Unfortunately, the effect of the 
light quark masses spoil this possibility [2]. While in lattice simulations it is in 
principle possible to reduce the light quark masses below their physical values and 
get arbitrarily close to the second order critical point, in heavy ion experiments we 
must live with a pion which has a mass comparable to the transition temperature. 
Near T^, the correlation length in the pion channel is shorter than T^^ [2], and as 
a result the misaligned regions almost certainly do not contain sufficient energy to 
radiate large numbers of pions. 

Here, we consider an idealization which is in some ways opposite to that of 
thermal equilibrium, that is the occurence of a sudden quench from high to low 
temperatures, in which the (o", tt) fields are suddenly removed from contact with 
a high temperature heat bath and subsequently evolve mechanically according to 
zero temperature equations of motion. In a real heavy ion collision, the phase 
transition proceeds by a process in between an equilibrium phase transition in 
which the temperature decreases arbitrarily slowly and a quench in which thermal 
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fluctuation ceases instanteously. Our goal is not to quantitatively model a realistic 
heavy ion collision as a quench. Rather, in studying the dynamics of the long 
wavelength modes of the chiral order parameter in a quench, our hope is that the 
qualitative behaviour in this model is representative of the physics which occurs 
in real processes in which the QCD plasma cools rapidly and is far from thermal 
equilibrium. 

We use the linear sigma model to describe the low energy interactions of pions: 
jC^ j d^^xf^^d^^rdi^cj^a - j{r<Pa-vY + Ha^ , (2.1) 

where A, v, and H oc ruq are to be thought of as parameters in the low energy 
effective theory obtained after integrating out heavy degrees of freedom. We treat 
(2.1) as it stands as a classical field theory, since the phenomenon we are attempting 
to describe is basically classical and because as a practical matter it would be 
prohibitively difficult to do better. 

Our numerical simulations of quenching in the linear sigma model are described 
in more detail in [1]. As initial conditions, we choose (f) and randomly indepen- 
dently on each site of a cubic lattice. Therefore, the lattice spacing a represents 
the correlation length in the disordered initial state. In [1] we made a crude at- 
tempt to choose initial distributions for and appropriate for a quench from an 
initial temperature of T = 1.2Tc. With initial conditions chosen, we then model 
the T — evolution of the system after the quench by evolving the initial config- 
uration using a standard finite difference, staggered leapfrog scheme according to 
the equations of motion obtained by varying (2.1). After each two time steps, we 
compute the spatial fourier transform of the configuration and from that obtain 
the angular averaged power spectrum. 

The central result of our simulations is that the power in the long wavelength 
modes of the pion field grows dramatically. While the initial power spectrum is 
white and while at late times the system is approaching a configuration in which the 
energy is partitioned equally among modes, at intermediate times of order several 
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times m~ the low momentum pion modes are oscillating with large amplitudes. 
When we used gaussian initial distributions for and with width v/2 and v 
respectively, the power in modes with k = 0.2a"^ ~ 0.3m7r is more than 1000 
times that in the initial white power spectrum. For initial conditions chosen to 
model an initial temperature of 1.2Tc, the amplification is less, but is still of order 
100. 

The amplification of low momentum modes which we observe in the numerical 
simulations can be qualitatively understood by approximating (f)^(f)a{x,t) in the 
equations of motion by its spatial average. After doing the spatial fourier transform, 
the equation of motion for the pion field becomes 

^n{k, t) = -{- A^;2 + A(02) (t) + e}7f{k, t) (2.2) 

where {4>^){t) means simply the spatial average of cfP'. At late times, {(jP") becomes 
time independent and takes its vacuum value, and the quantity in brace brackets 
in (2.2) becomes simply + A;^. Immediately after the quench, however, (0^) 
varies with time, and there are periods when (0^) < v^. A wave vector k mode of 
the pion field is unstable and grows exponentially whenever (0^) < — fc^/A. As 
(0^) varies, longer wavelength modes are unstable for more and for longer intervals 
of time, and, in agreement with the numerical simulations, are amplified relative 
to shorter wavelength modes. 



3. Charge separation does not occur 

The striking prediction (1.2) for the probability distribution of the neutral pion 
fraction R naturally leads to the question of whether there are similarly unusual 
fiuctuations in the electric charge itself, i.e. in the ratio of 7r+ to tt^ mesons. 
Formulae similar to (1.2) hold for the real fields t:^ — -^{t^^ + and tt^ = 
T^(7r"'" — 7r~), but not for tt"*" and 7r~. While the total electric charge must be 
conserved, there is no conservation law prohibiting the separation of charge into 
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regions of net positive and negative charge. We must determine whether long 
wavelength oscillations of the electric charge density grow. The charge operator 
Jo = 7i"^J^vr^ — TT^^TT^ measures rotary motion in the tt^ — tt^ plane. However, 
the amplification mechanism which operates following a quench kicks tt radially 
outward, and does not induce rotary motion. This heuristic argument is borne out 
in the simulations. Long wavelength oscillations of the electric charge density are 
not amplified. 

Notice that the question of whether charge separation occurs is a dynamical 
one, and has a straightforward dynamical answer for the quench mechanism of 
generating regions of misaligned condensate. In earlier work, Kowalski and Taylor 
[3] imposed isospin symmetry by hand in order to avoid the possibility of charge 
separation, which they consider physically implausible. The coherent states we 
reach are not isospin singlets, and we see no reason to impose that condition; 
nevertheless the intuition of Kowalski and Taylor is vindicated and there is no 
charge separation. 

4. How large are the regions of coherent pion oscillations? 

This question is of crucial phenomenological interest. In order to be observable, 
these regions must evolve into sufficiently many pions. We must ask, therefore, 
what are the longest wavelength modes of the pion field which get amplified? 
In our simulations, the answer is unequivocal — the wavelength of modes which 
are amplified is limited only by the size of the lattice on which the simulation is 
run. Alas, it is much harder to determine what will happen in a real heavy ion 
collision. The one thing that can be said with certainty is that the effect is not 
cut off by the inverse pion mass. Modes with k < tjIt^ are amplified. This is 
in marked contrast to the situation obtained in thermal equilibrium, and is why 
the phenomenon will only be detected if the plasma in a heavy ion collision is far 
enough from thermal equilibrium that quenching is an appropriate idealization. 
If is not the long wavelength cut-off, what is? The most optimistic (and 



6 



perhaps implausible) possibility, which is suggested by a literal interpretation of 
our simulations, is that coherent oscillations of the pion field in regions as large 
as the transverse extent of the plasma are possible. (The rapid expansion in the 
longitudinal direction will damp the growth of modes with k parallel to the beam 

— * 

relative to those with k transverse.) In a real collision, the size of modes which 
grow could perhaps be limited by the time available before the pions no longer 
interact and therefore can no longer be described by oscillations of a classical field, 
or perhaps by the size of regions of the plasma in which the energy density is 
reasonably homogeneous. At this point, the most that can be said is that long 
wavelength pion oscillations are amplified after a quench, and that their size is not 
limited by any microphysical length like but is limited dynamically, perhaps 
only by the system size. 



5. Outlook 



Although we have made many idealizations and approximations, it seems pos- 
sible that the essential qualitative feature of the phenomenon we have elucidated 
— long wavelength pion modes experiencing periods of negative mass^ and conse- 
quent growth following a quench — could occur in real heavy ion collisions. Given 
the explicit symmetry breaking which gives mass to the pions, one might have 
expected the dynamics following a quench to be featureless. The mechanism here 
discussed provides a robust counterexample. We have not come close to modelling 
a real heavy ion collision. While our treatment can surely be improved, it seems 
doubtful that quantitative theoretical predictions for heavy ion collisions will be 
possible. At the end of the day, the question of whether or not long wavelength 
pion oscillations occur will be answered experimentally. If a heavy ion collision is 
energetic enough that there is a central rapidity region of high energy density and 
low baryon number, and if such a region cools rapidly enough that the process can 
be modelled as a quench, this will be detected by observing clusters of pions of 
similar rapidity in which the neutral pion fraction R is fixed. This ratio will be 
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different in different clusters and will follow a distribution like (1.2). Were such 
a signature to be observed experimentally, it would be clear evidence for an out 
of equilibrium transition from a QCD plasma in which the chiral order parameter 
was initially disordered. 

REFERENCES 

1. K. Rajagopal and F. Wilczek, to appear in Nucl. Phys. B, |hep-ph/9303281| , 
Princeton preprint PUPT-1389, IASSNS-HEP-93/16, 1993. 

2. K. Rajagopal and F. Wilczek, Nucl. Phys. B399 (1993) 395. 

3. A. Anselm and M. Ryskin, Phys. Letters B226 (1991) 482; J.-P. Blaizot 
and A. Krzywicki, Phys. Rev. D46 (1992) 246; J. D. Bjorken, Int. J. Mod. 
Phys. AT (1992) 4189; J. D. Bjorken, Acta Phys. Pol. B23 (1992) 561; K. 
L. Kowalski and C. C. Taylor, preprint hepph/9211282, 1992; J. D. Bjorken, 
K. L. Kowalski, and C. C. Taylor, SLAC preprint SLAC-PUB-6109, 1993. 



8 



